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ABSTRACT 

The X-ray observation of AM Her in a very low state was performed with Suzaku in October 
2008. One flare event with a time scale of ~ 3700 sec was detected at the X-ray luminosity 
of 6.0 x 10 29 erg sec -1 in the 0.5 - 10 keV band assuming at a distance of 91 pc. The X-ray 
spectrum is represented by a thermal plasma emission model with a temperature of 8.67j^'^4 
keV. During the quiescence out of the flare interval, Suzaku also detected significant X-rays at 
a luminosity of 1.7 x 10 29 erg sec -1 in the 0.5 - 10 keV band, showing a clear spin modulation 
at a period of 0.1289273(2) days at BJD 2454771.581. The X-ray spectra in the quiescence 
were represented by a MEKAL + Power Law (PL) model or a single CEMEKL model, which 
are also supported by phase-resolved analyses. A correlation between the temperature and the 
volume emission measure was found together with historical X-ray measurements of AM Her in 
various states. In order to account for a possible non-thermal emission from AM Her, particle 
acceleration mechanisms in the AM Her system are also discussed, including a new proposal of a 
shock acceleration process on the top of the accretion column. 



Subject headings: novae, cataclysmic variables - stars: individual (AM Herculis) 
of particles 



plasmas - acceleration 



1. Introduction 

Since white dwarfs (WDs) exist everywhere in 
our Galaxy and have so large number density 
that one third of all the ast ronomical galactic 
objects are WDs (j Alien! 1 19731 ). they have some 
possibilities to make important contributions to 
unresolved long-standing mysteries in the astro- 
physics, such as the origin of Cosmic-ray parti- 
cles, Galactic ridge emissions, and so on. There- 
fore, it is important to understand the nature of 
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WDs in detail. Magnetic cataclysmic variables 
(MCVs) are binary systems which consist of late 
type stars and WDs with strong magnetic field 
strength of over 10 5 G. Cold gases from a com- 
panion star via a Roche-Lobe overflow accrete 
onto the magnetic pole(s) of the WD, and re- 
lease their gravitational potentials at the standing 
shock on the WD surface. The temperature of the 
shock heated plasma, called an accretion column, 
reaches about 10 8 K, and the plasma emits hard 
X- rays via bremsst rahlung process (see the review 
bv lPatterso3l994( l. Polars are sub-class of MCVs, 
which have so strong magnetic field strength of 
> 10 7 G that spin periods of WDs are synchro- 
nized to the orbital periods. 

In this paper, we report on the X-ray ob- 
servation of the pr oto-type star of pola rs, AM 



Her, with Suzaku (jMitsuda et al.l 120071 1. The 
distance to the object d, the magnetic field 
strength on the surface of the WD, Bwd, the 
mass of the WD, Af-yvD, its radius, -Rwd, and 
the spin period, P S pin, are measured as d = 91 
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c (iGansicke et all Il995l) B wr> = 13-30 MG 



(Schmid t et al.lll98U 11 983: Wickra masinghe et al 
Il984l : IWickramasinghe fe Martini Il985l). M W p 



0.39- 1.2 May (lYoung et al.lll98ll:lMukai and Chariesl 



1987t lBailev et al.l|l988HCropper et al.lll9 981. and 



3.09 hr IjStockman et al.l Il977l ). respec 



± spin 

tively, where Mq and Rq are the mass and ra- 
dius of the Sun, respectively. Due to changes of 
an accretion rate, M, from the companion star, 
the optical flux of AM Her shows two states; 
the high state with V magnitude of 13.5 mag 
and the low state with 15.5 mag. The previous 
X-ray observa tions in high states ar e performed 
by HEAO-1 (iRothschild et all Il98lh. EX OSAT 
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Osborne et al.l 1 19861) Ginaa (llshida et al.l [l991 



ASCA (llshida et al.l Il997i: iTerada et al.l (20041) . 
They all report the detail measurements of the 
thermal emission from the accretion column of 
AM Her. On the other hand, X-ray observations 
i n the low state a re reported only by BeppoSAX 
(|Matt et al.ll2000h . and those in an intermediate 
state are covered by Chandra ( Girish et al.ll2007h . 

X-ray studies of thermal emissions in a low 
state can help us to understand the physics of 
the accretion flow by exploring one of the impor- 
tant physical parameters, M , in a wide range. In 
addition, in a very low state, it should be easier 
to search for unknown faint emissions behind the 
thermal radiation. One of possible new-type com- 
ponents is a non-ther mal emission as suggested by 
Terada et al. (2008b) from Suzaku observations of 



Beardmore et alj|1995t ICropper et al.lll998l). an d 



an intermediate polar AE Aqr, which belongs to 
a sub-class of MCVs having weaker magnetic field 
strength than polars. Although the observation 
was actually not in a low state, the high-sensitive 
wide-band spectroscopy and phase-resolved analy- 
ses enabled them to find spiky pulsations like neu- 
tron star (NS) pulsars in the spin profiles with two 
independent instruments on-board Suzaku. They 
claimed that these spikes may have non-thermal 
origin from detailed analyses of X-ray spectra. 
Thus, the object is no w called as a "NS pul- 
sar equivalent of WD" (|de Jagerlll994l ) or just a 
"WD pulsar". X-ray deep observations in very 
low states could provide such a discovery; a sec- 
ond "WD pulsar" , which is important in studies of 
searches for a new cosmic-ray acceleration system. 
In this paper, we describe the X-ray observation 
of AM Her in a low state in section [21 summarize 



the results of Suzaku analyses during the flare and 
the quiescence in section [3j and discuss in section 

SI 

2. Observation and Data Reduction 
2.1. Suzaku observation of AM Her 

Suzaku is the fifth series of the Japanese X- 
ray satellite, carrying two X-ray in struments; the 
X-ray Imaging Spectrometer (XIS; iKovama et al 



| 2007h and the Hard X-ray Detector (HXD; lTakahashi et al 
120071 ). It has good capabilities to provide very low- 
background observations with high sensitivities in 
the 0.2 - 600 keV bandpass. Thus, it is suitable 
to search for a hidden emission from a faint object 
with a wide-band spectroscopy. In addition, high 
timing capability with the time resolution of 8 sec 
and 64 ysec for the XIS and the HXD, respec- 
tively ( Terada et al. 2008al ). can help us to follow 
the fast variability of the object. 

We observed AM Her with Suzaku from 2008 
October 29 20:22 UT to November 1 07:18 UT on 
the HXD-nominal position (OBSID=403007010; 
hereafter AMHER). According to optical reports 
by the American Association of Variable Star Ob- 
servers (AAVSO) 0, this observation date cor- 
responds to the low state of the object lasting 
from September 2008 to June 2009. In addi- 
tion to the main observation tagged as AMHER, 
we also performed an off-axis observation near 
the object to have the best use of the Suzaku 
sensitivities from 2008 November 1 07:18 UT to 
November 2 08:15 UT (OBSID = 403008010; here- 
after AMHER_OFF_AXIS). The off-axis obser- 
vation was carried out towards the position of 
(RA, DEC)=(282.08, 47.99), where no hard X-ray 
sources are contaminated in the field of view of 
the HXD PIN and GSO. On both occasions, the 
XIS was operated in the normal clocking mode 
without window/burst options bu t with the Space 



Char ge- Injection (SCI) function (jNakaiima et al 



2008). The HXD was operated in the nominal 



mode; half of 64 PIN diodes were operated with a 
voltage of 400 V and the others with 500 V. 

2.2. Data Reduction 

We used the data sets by the standard Suzaku 
pipe-line processing version 2.2.11.22, with the cal- 
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Table 1 

Source list in the background region. 



OBSID 



RA 



DEC 



Flux* 



index* 



Identification 



403007010 18: 

403007010 18: 

403008010 18: 

403008010 18: 



10.379 
34.923 
04.720 
17.747 



-49:47:06.95 
-49:48:32.84 
-47:57:34.24 
-47:52:25.25 



4-80l}f 3 
1 90+ ' 49 

l.»U_g 4 g 

°-'°-1.13 



2.41 



+ 1.27 
-2.40 



-, OO+0.34 
1 -°°-0.31 
o C7+0.67 
Z '°'-0.65 
0.31 
0.32 
+ 1.42 
1.26 



2.07 
2.44 



1WGA J1817.1+4947 
SUZAKU J1816.6+4948.5 § 
SUZAKU J1848.0+4757.6 § 
SUZAKU J1848.3+4752.4§ 



f X-ray flux in 0.5 - 10 keV band, 10~ 14 erg sec -1 cm 2 . 
X Photon index of the object. 
§ Newly discovered with Suzaku in this work. 



ibration version (CALDBVER) of hxd20081009, 
xis20081009, xrt20080709, and xrs20060410. We 
used the ftools in HEADAS 6.6 with XSPEC ver- 
sion 11.3.2ag. 

The source was detected with the XIS at the 
average count rate of 0.03 - 0.04 cnt sec -1 per 
sensor in the 0.5 - 10 keV band. Cleaned events 
of the XIS data were obtained with the stan- 
dard criteria of the pipe-line process. The ex- 
posures of AMHER and AMHER.OFF.AXIS ob- 
servations were 108.5 ksec and 44.3 ksec, respec- 
tively. On-source events of the XIS were then ob- 
tained by accumulating within 6 mm (4.3') ra- 
dius from the image centroid of the object in 
the AMHER dataset. Several faint sources (in- 
cluding one known object) were detected both in 
AMHER and AMHER_OFF_AXIS datasets with 
the XIS count rate of less than about 0.005 cnt 
sec -1 , which is comparable to the count rate of 
the background and corresponds to the X-ray flux 
of 0.7 x 10 -13 erg sec -1 cm 2 as listed in ta- 
ble [1] Thus, the background events were accu- 
mulated by source free regions of AMHER and 
AMHER_OFF_AXIS, where > 6 mm outer than 
these faint sources or AM Her. We summed up 
these two background spectra, because these spec- 
tral shapes and the flux of the two blank sky 
datasets were consistent within 5 % in the 0.5 - 
10 keV range. 

In the data reduction of the HXD, we first 
reprocessed the unscree ned events by the tool 
'hxdpi' and 'hxdgrade' (jTerada et al. 2005 ) us- 
ing the CALDB files (including the GSO gain- 
history file named ae_hxd_gsoght_20090131.fits to 



correct the PMT gains of GSO), and then ob- 
tained cleaned events of PIN and GSO with the 
standard criteria of the process. The exposures 
of the HXD of AMHER and AMHER.OFF.AXIS 
were 94.4 ksec and 40.3 ksec, respectively. The 
non X-ray background (NXB) events of PIN and 
GSO are estimated and provided by the HXD 
team with the methods by Fukazawa et al. (2009). 
We used NXB events of both PIN and GSO with 
METHOD='LCFITDT (bgd_d)' and version of 
METHODV='2.0ver0804'. The systematic errors 
of the NXB can be checked by comparison be- 
tween the actual data and the NXB data dur- 
ing the night earth observations. In AMHER 
and AMHER_OFF_AXIS observations, from night 
earth spectra with exposures of 38 ksec and 59 ksec 
for PIN and GSO, respectively, the reproducibili- 
ties of PIN and GSO NXBs were about 4 and 1 %, 
respectively. Within these systematic errors, the 
events detected with PIN were consistent with the 
Cosmic X-ray background (CXB) level, and those 
with GSO were consistent with being no signals. 

3. Analyses and Results 
3.1. X-ray Light Curves 

Fig.[T]shows background-subtracted light curves 
of AM Her, obtained with the XIS and the HXD. 
The arrival times in this figure are barycen- 
tric corrected with the tool named 'aebarycen' 
2008al ). In the XIS data, a flare 



(Terada et al 



phenomenon starting from BJD 2454770.78 was 
observed with a duration of about 7 ksec. There 
was not any brightening in the background data 
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of the XIS, and so the flare should arise from 
AM Her itself. On the other hand, the PIN data 
showed no significant flare synchronized to the XIS 
one. Here, we defined the duration between BJD 
2454770.777416 - 2454770.863308 as an epoch of 
the flare, and the others as the quiescence. The 
effective exposures of the XIS during the flare and 
the quiescence were 5.0 ksec and 103 ksec, respec- 
tively. 
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Fig. 1. — Background subtracted light curves 
taken with the Suzaku XIS and HXD in 0.5 - 1.5, 
1.5 - 4.0, 4.0 - 10.0, 13 - 40 keV band, from top 
to bottom, respectively. 
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Fig. 2. — Light curves folded at the spin period of 
Pxis = 0.1289273 days shown in 0.5 - 1.0, 1.0 - 
2.0, 2.0 - 4.0, 4.0 - 7.0, 7.0 - 10.0, and 13 - 40 keV 
band from the top to bottom panels, respectively. 
The phase 0.0 corresponds to BJD 2454771.581. 
The upper six panels are taken with Suzaku, and 
the lowest panel sh ows the same plot but with 
ASCA CIS in 1993 dTerada et alJl2004h . 



We search for the spin period from the XIS data 
in the 0.5 - 10 keV band during the quiescence, 
and obtained the period P X is = 0.1289273(2) 
days at the epoch BJD 2454771.581. The value 
is con sistent with the P S nin b y previous qbserva - 



tions (|Stockman et al.l Il977t Kafka et all |2005) 



Ps 



xis = Pspin- Here, the epoch was simply 
determined by the peak time of the pulse pro- 
file fitted with the sine function. Fig. [2] shows 
the energy-resolved light curves folded at Pxis- 
Clear spin modulations are observed especially in 
the soft energy band below 4 keV. They have 
similar pulse s hapes to those taken with ASCA 
(jTanaka et al.l 119941 ) during the high state, al- 
though significant energy dependence of the spin 
amplitudes is detec ted, unlike those in high states 
(jlshida et al.lll997l ). 



3.2. X-ray Properties during the Flare 

Fig. |3] shows the light curve of the XIS in the 
0.5 - 10 keV band around the flare. The folded 



0.5 - 10 keV light-curve in the quiescence (like Fig. 
[2]) are overlaid periodically for comparison. If wc 
subtract this spin variation, X-ray counts around 
the phase bottom (BJD 2454770.86 - 2454770.88) 
became almost zero, but the others showed an ex- 
ponential decay. Then, we fitted the light curve 
without the phase-bottom epoch by the exponen- 
tial function of exp((T — To)/7"flare), where T, To, 
and Tfl are are time, starting time of the flare, and a 
decay constant of time, respectively. We obtained 
the time origin of T = BJD 2454770.76 ± 0.03 and 
the decay constant of Tfl arc = (4.3 ± 0.2) x 10~ 2 
days. 

The X-ray spectra during the flare are shown 
in Fig. |U Since there were no significant de- 
tections by the PIN and GSO within system- 
atic errors of the NXB (section 12. 2p , only XIS 
events were plotted in the figure. First, we tried 
to reproduce the X-ray spectra by an emission 
model from an optically-thin thermal plasma, 
which can be generated at the post-shock region 
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Table 2 

Best fit model parameters in the flare epoch of AM Her. 



Model 


n H 
10 20 cmr 2 


kT 
keV 


a * abundance 
solar 


Fluxt 
10~ 12 erg cm~ 2 s _1 


xl (d.o.f) 


phabs * MEKAL 
phabs * CEMEKL 


< 1.3 

< 4.0 


17 6+ 13 - 9 


u - '°-0.26 

l 64+ - 70 o qn +0 - 40 

^ -0.65 u - yu -0.35 


D ' UO -0.38 


1.00 (205) 
0.99 (203) 



f X-ray flux in 0.5 - 10 keV band. 
\ The power of DEM, as present ed by DEM oc (T/T m a X ) a " 1 rf(T) in the CEMEKL model 

dDone fc Osbornlll997h . 



of the accretion gas-flow onto the WD via en- 
ergy releases of their gravitational potential. The 
spectra obtained with Suzaku were well repro- 
duced by s uch a thin-therma l plasma model called 
MEKAL dMewe et al.l [l985t bedahl et all Il995t 
iKaastra et al.l 1996 ). as shown in tableland Fig. 
H The fitting to the XIS-0, 1, and 3 datasets 
were performed simultaneously, and cross nor- 
malizations were consistent with the value of 1.0. 



Although the temperature, kT = 8.67 



+ 1.31 

-1.14 



keV, 



obtained during this very low state, was slightly 
lower than t hose in high states (e.g., 13.5 keV 
or 11 keV bv lBeardmore et ah 1995 ; Ishida et~al"1 



1997L respectively) the metal abundance, 0.761°' Jg 



solars, was consistent with values in high states. 
Thus, the hot plasma would come from the ac- 
cretion matters from the companion star. The 
spectra requires no significant photo-electric ab- 
sorption as shown in table [2] this fact is also 
self-consistent with the situation of a very low 
accretion rate. 

You may notice that there remains a line struc- 
ture around 2.0 keV in the residual plot of Fig. [4j 
although they were insignificant. If it was real, 
the features should come from a light element 
Si, from which emission should be suppressed in 
such a high temperature as kT = 8.67 keV. For 
a reference, we further tried an one-dime nsional 
coolin g-flo w mod e l as or iginally proposed bv lHoshj 
(|l973h and|Ajzu| (|l973h . and found that no more 
free parameters improve the fitting statistically. 
The spectra were also well fitted by such a model 
as the CEMEKL model jDone & Osbornl Il997t 



Baskill et al.l 2005), but there was no large im- 



provement in the chi-squared value, as shown in 
table [21 In this fitting, the shock temperature, 



T max , were obtained as 17.6+g 3 6 9 keV, which is 
higher than the average temperature by the sin- 
gle MEKAL fitting, because T max represents the 
maximum temperature just below the shock front 
in the accretion plasma. The value itself was rea- 
sonable compared with the gravitational potential 
of the WD. The power of the temperature gradi- 
ent, a, of the differential emission measure, DEM, 
was a = 1.64+ ,'!?, which was different from the 



theoretical expectation of 0.5 (jlshida et al.l 11994') 



with a simple assum ption of an one-dimensional 
stable cooling-flow bv lAizul (Il973h 



3.3. Average Spectral Analysis in the Qui- 
escence 

The X-ray spectra during the quiescence with 
Suzaku were shown in Fig.[5j Cross normalizations 
between XIS-0, XIS-1, and XIS-3 were consistent 
with 1.0, as already checked in section 13.21 and 
thus datasets of all the XIS chips were summed in 
the figure. As implied in section l3~T1 and Fig. [21 the 
spin pulse profile during the quiescence shows the 
similar modulations as in the high states, soft X- 
ray emissions should have thermal origin. There- 
fore, we tried to represent the spectrum by a single 
MEKAL model, following the steps on the analy- 
ses during the flare (section I3.2[) , but the fitting 
was not acceptable (x 2 = 1.30) as shown in table 
[U The metal abundance was quite low as 0.06+Q Q2 
solar, and so the best fit model has no physical 
meanings. 

In the hard X-ray band, the spin pulse profiles 
(Fig. [2]) showed shallower modulations than those 
in the softer energy band; the pulse shape was 
consistent with being flat in above 10 keV range. 
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Fig. 5. — X-ray spectra with Suzaku during the quiescence of AM Her in the very low state. The crosses 
represents the XIS data summed of XIS-0, XIS-1, and XIS-3. Lines show the best fit models of single 
MEKAL, double MEKAL, MEKAL + power law, and CEMEKL, in top left, top right, bottom left and 
bottom right panels, respectively (see the text). Parameters of the models are listed in table G3 



Thus, in the spectral analyses, we then added an- 
other component in the hard X-ray band, repre- 
senting another thermal emission (case-1) or non- 
thermal component (case-2). 

In the first case that hard X-rays have a thermal 
origin (case-1), we tried a double MEKAL model 
to represent the spectrum. As a result, the fit- 
ting became acceptable (xl — 0.76), as shown 
in table [3] and Fig. [5j However the tempera- 
ture of the higher component was quite high as 
> 17.7 keV and the metal abundance was still too 



low as 0.061 



0.02 



solar, compared with the result 



during the flare (section 13. 2[) and /o r reports by 
previous observations in high states (llshida et al 



19971 : iTerada et all 12004 iGirish et all 120071 ). In 
the fitting, the two parameters of the abundance 
in both MEKAL components were linked, but the 
situation did not change even when we set these 
two parameters free from each other; 1.03lQ Qg 
keV plasma with abundance of 0.06to'oi solar and 
> 14.8 keV plasma with < 0.85 solar abundance 

(xl = i.oi). 

In the second case that hard X-rays have a 



non-thermal origin (case-2), we tried MEKAL 
+ PL model and obtained an acceptable result 
(xl = 0.61), as shown in table [3] and Fig. [3] 
The metal abundance of the MEKAL component, 
0.80^15 solar, became consistent with previ- 
ous observations (Ishida et al. 1997: Terada et al 



20041 : iGirish et al. 2007) and this work during the 



flare (section 13.21) . The photon index of the non- 
thermal component was 2.36^014- The X-ray flux 
in the 0.5 - 10 keV band of the non-thermal com- 
ponent was (1.5 ± 0.1) x 10~ 13 erg sec -1 cm~ 2 , 
which corresponds to the luminosity in the same 

energy band of (1.5±0.1) x 10 29 (g^j erg sec -1 . 
Note that a hump structure around 6-7 keV 
in the residual panel of Fig. [5] bottom left was 
statistically insignificant; if we tentatively add a 
Gaussian model here to MEKAL + PL model, 
we obtained the central energy of the Gaussian 
at 6.3^g'g keV, where the improvement of xl was 
only 0.02 corresponding to the F-test value and 
the probability of 1.76 and 0.18, respectively. 

As an application of case-1 fitting (double 
MEKAL model), we also tried the multi-color 
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MEKAL+ Power Law 



MEKAL + MEKAL 



MEKAL + Power Law 



CEMEKL 
---it, 




MEKAL + Power Law 



MEKAL + MEKAL 



MEKAL + Power Law 



Fig. 6. — Left and right panels show the X-ray spectra in the top and bottom phases, respectively. Crosses are 
the XIS data. The top panel in each figure shows the spectrum with the best fit model of MEKAL + power 
law, and bottom panels are the residuals of data from models, MEKAL, MEKAL+MEKAL, MEKAL+Power 
Law, and CEMEKL, towards lower panels. Parameters of the models are summarized in table [3J 



plasma model, CEMEKL, as already used in sec- 
tion 13.21 (case-3), and got an acceptable result 
(xl = 0.32) as summarized in table [3] and Fig. 
[5j The metal abundance and the fcT mav of the 
best fit model were reasonable as 0.48 



+0.22 
-0.21 



so- 
lar and 6.37l*'3y keV, respectively. The power 
a of the Volume Emission Measure (VEM) was 
a = 0.58^0 ^, which was c onsistent with the t he- 
oretical value of a — 0.5 l lshida et al.lll994T) in 
an ideal situation of the one-dimensional steady 
cooling flow. In summary, the phase-averaged X- 
ray spectrum of AM Her in the quiescence were 
reproduced not by single MEKAL model, but well 
by three kinds of models; (case-1) MEKAL + 
MEKAL, (case-2) MEKAL + PL, and (case-3) 
CEMEKL, although case-1 requires an unrealisti- 
cally low metal abundance. 

3.4. Phase resolved analyses in Quiescence 

In order to restrict the emission model in the 
quiescence, by a combination of results of the tim- 



ing analyses (section 13.11) and the spectral ones 
(section I3.3[) . we compared the best fit models 
of the phase-resolved spectra with the spin pulse 
modulation in Fig. [21 quantitatively. 

First, we tried phase-resolved spectral analyses. 
From the spin profile in Fig. [21 we defined the top 
and bottom phases as epochs at spin phases of cf> — 
0.59375 - 1.21875 and 0.21875 - 0.59375, respec- 
tively, and selected XIS events in these two phases. 
The X-ray spectra in both phases are shown in Fig. 
El Following the phase averaged studies in section 
[331 we tried to fit these spectra with (0) MEKAL, 
(1) MEKAL+MEKAL, (2) MEKAL+PL, and (3) 
CEMEKL models, and got results shown in table 
[3] and Fig. As a result, there were no strict 
differences statistically among these models ex- 
cept for case-0, single MEKAL model, as already 
shown in the phase averaged analyses. In addi- 
tion, the X-ray spectra of the pulsation component 
(i.e., subtraction from the spectra during the top 
phase of those during the bottom phase) were also 
well represented both by thermal and non-thermal 
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Table 3 

Best fit model parameters during the quiescence of AM Her. 



Model 


fcTi 


kT 2 or fcT max or F* abundance Flux § 


xl (d.o.f) 




keV 


keV solar 





Phase Average 




















MEKAL 


i 90+0.09 
1 - zo -0.13 






n nfi+ 003 


1 


1 9+0.08 
lo -0.11 


1 


.30 


(51) 


MEKAL + MEKAL 


1 02+ - 07 


> 11.1 




0.06±0.02 


1 


7O+0.21 
■ ' y -0.30 





.76 


(49) 


MEKAL + PL 


1.10±0.09 




2-36±g;« 


o.8ol8:i8 


1 


79+1-27 
,o -0.08 





.61 


(49) 


CEMEKL 




c r 7 +3.63 
D - ' -2.57 


o.58i8:S 


0.4818:2? 


1 


7fi +0.20 
■ '"-0.54 





.32 


(50) 


Top Phase 




















MEKAL 


1.0718°? 






0.05±0.02 


1 


C9+0.14 
uz -0.13 


1 


.24 


(51) 


MEKAL + MEKAL 




> 6.2 




0.04±0.02 


2 


9R+0.14 

oo -0.48 





.73 


(49) 


MEKAL + PL 


1-0418:85 




2 42+ - 17 

z -^ z -0.18 


o.8ol8:| 


2 


19+1.85 
■ lo -0.08 





.67 


(49) 


CEMEKL 




6.37t?; 3 3 4 7 


U.OO_ 33 


o.5ol8: 5 2 I 


2 


-i 7+0.23 
1 ' -0.62 





.37 


(50) 


Bottom Phase 




















MEKAL 


l-68±C;I? 






u - lq: -0.10 





Q7+0.86 
y ' -0.22 





.24 


(51) 


MEKAL + MEKAL 


i 9.7+0.59 
i -°'-0.35 


79.9l° 79 1 8 




o.nl8:S 


1 


97+0.45 
■ z ' -0.31 





.20 


(49) 


MEKAL + PL 


1 - 4d -0.36 




9 /~i7+6.66 
z - u ' -1.65 


0.38181? 


1 


99+0.26 
■ zz -0.14 





.19 


(49) 


CEMEKL 




5.76t 2 3 3 ° 


1 - zy -3.50 


n 7R+ 317 
u./o_ 67 


1 


2fi +ln 

ZD -0.56 





.17 


(50) 




f Valid for 


CEMEKL model, the power of DEM. 













\ Valid for Power law model, the photon Index. 
§ The X-ray flux in 0.5 - 10 keV band in 10~ 13 erg cm~ 2 s _1 . 



Table 4 

Best fit model parameters of the pulsation spectrum during the quiescence. 



Model 


kT 
keV 


r* abundance 
solar 


Flux§ 


Xl (d.o.f) 


PL 

MEKAL 
CEMEKL 


0.78±0.17 

10 8 +89 2 

lu -°-10.8 


9 Q9+0.41 

o.o5l8:° 8 3 

0.121™ • • o.84i8:a 


i-oo!8:S 
o.73l8: 2 ! 

i-ool8: 78 


0.47 (52) 
0.32 (51) 
0.27 (50) 


f Power of the DEM in the CEMEKL model, a. 
\ Photon Index of the PL model. 
§ X-ray flux in the 0.5 - 10 keV band in 10~ 13 erg cm~ 2 s _1 . 
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AM Her, Flare 

amher_flare_fit02b, nH<1.27e20, kT=B.67-1. 14+1.31 keV, Ab=0.76-0.26+0.28 



2454770.5 2454770.7 2454770.9 2454771.1 
BJD 

Fig. 3. — Top panel shows the light curves of 
the XIS in the 0.5 - 10 keV band near the flare 
event shown in black crosses, and pulse profiles of 
the XIS in quiescence folded at the spin period 
of PxiS = 0.1289273 days are also shown in gray. 
The dashed lines show the average count of the 
quiescence and the best-fit exponential curve (see 
the text). Bottom panel represents the hardness 
ratio of the XIS; counts in 1.5 - 10.0 keV band 
divided by those in 0.5 - 1.5 keV band. 

models as shown in table 2] i.e., we could not dis- 
tinguish these origins statistically, only from the 
phase-resolved analyses. 

Next, we evaluated numerically the amplitudes 
of spin modulations as a function of the energy. 
Since the modulations are well reproduced by a 
single sine function (reduced x 2 /d.o.f ~ 0.6) and 
power of higher harmonics of the waveform arc 
negligible in estimating amplitudes, we fitted the 
energy-resolved pulse profiles in Fig. [5] by a sine 
function with a constant model, and derived the 
ratios of the pulsation amplitudes to the average 
ones (i.e., a DC component). As shown in the 
results in Fig. [7] crosses, it was clearly seen that 
there were shallower modulations in higher energy 
band. This situation is different from those seen 





■ 1 ■ ■ 1 







channel energy (keV) 

Fig. 4. — X-ray spectra with Suzaku during the 
flare of AM Her in the very low state. The black, 
red, and green crosses represents the data of XIS- 
0, XIS-1, and XIS- 3, respectively. Lines show the 
best fit model of MEKAL with parameters in table 

m 



suit with phase-resolved spectral models, we also 
plotted in the same figure the ratios of models be- 
tween the pulsation components (i.e. subtractions 
of best fit models in the bottom phase from those 
in the top) and the best fit models in the aver- 
age (section ESI . Therefore, MEKAL + MEKAL 
model (case-1) requires a large modulation in the 
hard X-ray band, which was inconsistent with the 
manner of the spin modulation. The best model 
was MEKAL + PL (case-2) in Fig. although 
CEMEKL model (case-3) could not be excluded. 
Note that the negative values of the ratio seen in 
the MEKAL+PL model is due to the difference of 
the best-fit values of the PL index and it is con- 
sistent with being above 20 keV. 

4. Discussions 

4.1. Plasma Parameters during the Flare 

A flare event was observed during the X-ray 
observation of AM Her in the very low state with 
Suzaku, as shown in section 13.21 Same kinds of 
flares were reported in other polars . such as U Z 



For d uring a low X-ray luminosity iPandel et al 



with ASCA in the high state (jlshida et al.lll997l) 
Then, in order to check the consistency of this re 



([20021 ). The X-ray flare of AM Her decayed with a 
time scale of Tfl are ~ 3700 sec, and the X-ray spec- 
tra was well reproduced by the thermal plasma 
model with the average temperature of kT = 8.7 
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MEKAL + MEKAL 




Amplitude of pulse profiles 



Energy (keV) 



Fig. 7. — Crosses represent the amplitudes of the 
energy-resolved pulse profiles in Fig. [2] to the av- 
erage counts (see the text). Lines show the ratio 
of models taken by phase-resolved spectral analy- 
ses in table |3] to the average models; MEKAL + 
MEKAL, MEKAL + PL, and CEMEKL (Also see 
the text). 

keV (table [5]) and the volume emission measure 

(VEM) of 3.12 x 10 53 (g^) cm" 3 . In this sec- 
tion, we discuss on a rough estimation of plasma 
parameters during the flare, i.e., the electron den- 
sity, n e , and the plasma size (or the shock height of 
the plasma, h) , by using three observational values 
of r flarc ,fcT, and VEM. 

First, from the observational parameter of Tfl are , 
we can estimate the value of h, if we assume the 
flare time scale (rfl arc ) is mainly determined by 
the cooling time scale of the plasma, r coo i. This 
assumption may be correct because the temper- 
ature seemed to be gradually decreased through 
the flare as seen in the bottom panel of Fig. [3] 
According to Aizu (1973), the shock height is nu- 
merically described as h = 0.605 • u s ■ t coo \, where 
u s is a shock velocity of the plasma, which is about 
10 9 cm s _1 on a typical gravitational potential of 
a WD in a strong shock case, as 



1.2 x 10 



<) 



AfwD 

0.88M S 



1/2 



Rwd 
10 9 cm 



-1/2 



stf) 



Thus, we get 



h = 



3.1 x 10 



XX ( ''"cool 



V3700 s/ 



WD 



1/2 



R 



WD 



10 9 cm 



-1/2 



cm{2) 



4 0.88M Q/ 

This value is quite large in the binary system of 
AM Her compared to the binary distance, a, which 
is calculated by the third law of Kepler as 

1/3 / p \ 2/3 



a = 9.8x10" 



WD 



M 2 



0.88M n 0.3M C . 



\PxisJ 



cm, 
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wher e M 2 the mass of the companion star (| Bailey et al 
19881) and P or b is the orbital period of the binary. 



In other words, this solution requires that the bi- 
nary space is fully filled with the plasma, and it 
is unphysical. Therefore, the cooling time scale 
should be shorter than the flare time scale (i.e., 
Tcooi <C Tfl arc ) and Tflare is just described by a time 
scale of the accretion onto the WD from the com- 
panion star. This statement is supported by the 
observational fact seen in Fig. [3] top; the X-ray 
flux became faint during the epoch corresponding 
to the dim phase in quiescence, and thus the flare 
phenomenon occurs near the magnetic pole on the 
WD surface. 

Second, from the observational parameter of 
kT, we can constrain values of n e and h, under 
the condition that t coo j is not limited by Tfi are as 
discussed above. Here, we use n c as the elec- 
tron density just below the shock front, and we 
assume that the maximum temperature obtained 
with the CEMEKL model, A;T max , represents the 
shock temperature, fcT sri , of the multi-color accre- 
tion colu mn; i.e., we assum e fcT sr| = fcT max . Ac- 
cording to lLamb fc Masters! (|1979l ) , the cooling of 
the plasma occurs via free-free emission and cy- 
clotron radiation, and the time scale of the cooling 
is described as 

T cool = / _ 2 _ 2 i I 4 J 

where eg and e C y C are the volume-emissivities av- 
eraged along the accretion column of the free-free 
and cyclotron emissions, respectively. The aver- 
aged volume-emissivity of free-free emission are 
quantitatively described as 



i x 7.3 x 10 8 



kT f 



sh 



17.6 keV 

2 



1/2 



h x 10 15 cm- 3 
erg sec 



- 1 cm" 3 , 



(5) 
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where A is a correction factor of the plasma 
structure; when we put z as a position from 
the bottom of the plasma, and the plasma 
has a vertical dependence of kT oc nj 1 oc 
(z/h) p , then A is described as ^zsp- The c y~ 
cl otron emissivity e c ~~v C is given i n equ ation (10) 



of lWu, Chanmugam, and Shavivl ([19941 ) which de- 
pends on h - 85 . Here, the value h was determined 
by T coo i in equation (j4j, which has a dependency 
with e C y C again. Thus, we need to solve the simul- 
taneous equations between eg, e C y C , h, and n e to 
get self consistent values. We performed numer- 
ical calculations to solve them to derive a rough 
estimation between n c and h. In the calculation, 
we assume the Aizu solution (1973) where the cy- 
clotron cooling was ignored; i.e., we assume the 
relation h ~ 0.6u s r coo i and set the plasma struc- 
ture at p = 2/5. We also ignored the effects of the 
gradient of the gravitational potential of the WD 
in the accretion plasma. The solutions, when the 
magnetic field st rength is a typical v alue of AM 
Her B = 20 MG ((Schmidt et al.lll98lh . are shown 
in Fig. HI 



10 s 



i0 7 



to 6 



"2r-s.aj).. 

-acsje/i 




kT 1 ' = 17.6keV 
VEM = 9.6 x I0 5 W 3 




/ 





He (cm~ J ) 



10' 



Fig. 8. — Thick lines show numerical solutions 
of h and n e , based on calculations on kT s h with 
Suzaku during the flare phase of AM Her in the 
very low state. Dashed lines show the limitations 
from VEM with assumption of a = 2, 4, 8, 16. See 
the text. 

Finally, we can use the third observational 
value, VEM, to limit the parameters h and n c . 
To evaluate the volume of the plasma, we assume 
that the accretion column has a cylinder shape 
with a radius of r, which is a times size of the 
shock height (i.e., 2r = ah). Dashed lines in Fig. 



[8] represent the limitation from VEM values to h 
and n in cases of a = 2, 4, 8, and 16. In the plot, 
we used the VEM value from a simple MEKAL 
fitting; since the model assumes uniform density 
and the plot is shown in the shock density, we mul- 
tiplied a correction factor of 1/(1 — 2p), which is 
4 when we assume the Aizu solution (1973). As a 
result, we can get a rough values of h and Ti c as a 
function of a; for example, n c ~ 5 x 10 15 cm" 3 and 
h ~ 10 7 cm when a = 4. Although many solutions 
are allowed in this estimation, quite small values 
of a are not realistic from Fig. [SJ i.e., the plasma 
has a coin-like shape rather than tall-cylinder one. 

4.2. Structure of the Thermal Plasma in 
the Quiescence 

According to the results of the data analyses 
in sections 13.31 and 13.41 the X-ray spectra of AM 
Her during the quiescence can be described either 
by a multi-color-plasma CEMEKL model with a 
temperature of kT max = 6.7 keV, or by a combina- 
tion of thermal and non-thermal model of MEKAL 
with kT — 1.1 keV plus the PL component. In 
this section, we focus on the nature of the ther- 
mal component, whereas non-thermal phenomena 
in AM Her will be discussed in the next section 
1431 

First, we estimated the plasma parameters, n c 
and h, during the quiescence by performing nu- 
merical calculations as already shown in section 
14.11 The VEMs during the quiescence were one 
or two order-of-magnitude lower than that in the 
flare; VEM = 2.5 x 10 52 cm 3 (CEMEKL) or 
1.4 x 10 51 cm 3 (MEKAL+PL). To obtain the value 
of fcT sh , we assumed kT sh = /cT max for CEMEKL 



model as in section |4~T1 We cannot estimate kT s h 
directly from kT of MEKAL model, so we sim- 
ulated PHA spectra from CEMEKL model with 
various fcT max on XSPEC and fitted with MEKAL 
models to obtain a rough relation between kT and 
fcT max . In the simulation, we used the XIS re- 
sponse on the HXD nominal position and fitted 
spectra in the hard X-ray band of the XIS, 4-10 
keV. Then, we got a A:T max - kT relation as 



log(fcT max ) - 1.231og(*T) + 0.14. 



(6) 



Thus, the value kT = 1.1 keV by MEKAL analysis 
corresponds to fcT s h ~ 1.6 keV of the multi-color 
column. Using these VEM and fcT s h values, we 
performed the same numerical calculations again 



11 



as in section I4.1[ and got plasma parameters as 
n c ~ a few xlO 14 cm -3 and h ~ 10 7 cm when 
a = A. Although there remains an uncertainty on 
the plasma shape a, the column height seems to 
be almost the same order of magnitude as that 
in flares, and the electron density got one or two 
order of magnitude lower than in flares. 

Second, in order to compare the Suzaku results 
with other past observations in the high, medium, 
and low states, we plotte d kT and VEM values 
by several X-ray satellites dBeardmore et aL 19951: 



2004; 



Girish et al 



Ishida et al.lll997t iMatt et al.ll200d : ITerada et al 



2007|) in Fig. M In this figure, 



the kT values obtained with MEKAL models are 
corrected into fcT s h equivalents by the kT-kT max 
relation of equation ([6]). Therefore, observation- 
ally, the plasma temperature gets lower as VEM 
gets smaller, and the Suzaku points in the faintest 
state also follow this trend. Hereafter, we call this 
relation as the VEM - kT trend. 

Finally, to interpret the VEM - kT trend in Fig. 
|9] quantitatively, we performed further numerical 
calculations demonstrated in section 14.11 In Fig. 
HO we plotted contours of VEM and kT with vari- 
ous conditions of a, n c , and h. The calculations in- 
cludes the free- free and cyclotron cooling processes 
with several assumptions described in section 14.11 
Therefore, the results of calculations suggest that 
the changes on the plasma parameters occur only 
on n e , keeping the column shape almost constant, 
even when the VEM varies by two order of mag- 
nitudes. As shown in Fig. [9j one realistic set of 
parameters is a = 4 and h = 10 7 cm, which were 
already presented as an one possible solution of 
the flare plasma in section 14.11 Other possibili- 
ties are larger a at lower h = 10 6 cm or smaller 
a at higher h — 10 s cm, but the latter case is 
less feasible because analytic solutions vanish at 
small a values, as already mentioned in section 
14.11 In other words, our calculation does not ex- 
clude other interpretations that both n e and the 
plasma shapes varies as VEM changes. 

If the plasma shape is constant in various states 
keeping the same shock height h, then the grav- 
itational potentials at the shock front are not 
changed by its accretion rate. The variation of kT 
as observed in Fig. [§] may be caused by possible 
changes of the energy-conversion efficiency rj at the 
shock, however, the adiabatic index of the accre- 
tion gas is required to be less than 1 to change r\ by 



> 

1 1.0 




to 



■'(flare) 




1051 1Q52 1053 10W 1055 

VEM (cm-3) 

Fig. 9. — Crosses with circle show scatter plot 
of kT and VEM, obtained by X-ray observations 
of AM Her iBeardmore et all (Il995l); llshida et aT 



(Il997l); IMatt et all (|2000h : ITerada et all (|2004h : 



Girish et all (|2007l) . Lines show relations be- 
tween kT and VEM by numerical calculations on 
the various cases of h = 10 6 ,10 7 ,10 8 cm, n e 
10''. 10 15 , 10 16 cm -3 , and a = 2,4,8. Conditions 
are shown in the figure. 

factor two or three; this situation is not realistic. 
Another possibility is an observational effect; X- 
ray observations detect the free- free emissions and 
thus are essentially sensitive to the dense bottom- 
part of the column where the plasma is already 
cooled down mainly by the cyclotron radiation. At 
lower accretion rate, the effect will be significant 
because the cyclotron cooling becomes dominant 
in the lower densities as discussed in section 14.11 

4.3. Particle Acceleration on the Accre- 
tion Shocks 

In this section, we assume that the hard X-ray 
emission detected with Suzaku has non-thermal 
origin and we concentrate on the particle accel- 
eration processes in AM Her. As demonstrated 
in sections 13.31 and 13.41 the X-ray luminosity of 
the power-law component of AM Her was L Xj pl = 

1.5 x 10 29 erg sec" 1 at d = 91 pc in 0.5 - 10 keV 
band, which is the same order of magnitude as 
that of the first WD-pulsar candidat e, AE Aqr, at 

6.6 x 10 29 erg sec -1 at d = 102 pc (|Terada et al ' 
2008bT ). Thus, neither bremsstrahlung nor curva- 
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ture radiations are feasible as a physical process 
of non-thermal emissions from a viewpoint of ra- 
diation efficiencies, although another possibility is 
the synchrotron radiation from about MeV elec- 
trons under the stro ng magnetic field of the WD 
(|Terada et al.ll2008bl ). 

If particles are accelerated by an electric po- 
tential induced by the spin rotation of the mag- 
netized WD up to 10 14 volts, the acceleration 
process should occur on the same place as the 
non-thermal emissions; i.e., near the WD surface 
where the magnetic field is strong enough to cre- 
ate hard X-rays via synchrotron process. How- 
ever, in such a polar-cap case, the emissions should 
present spin pulsations like AE Aqr case, whereas 
no modulations were observed in the hard X-ray 
band from AM Her (section 13. ip . Thus, elec- 
tric acceleration mechanism may be less feasible, 
although we need more theoretical studies. In 
this section, we further consider another pos sibil- 
ity of diffusive shock-acceleration m echanism (jBelll 
19781: iBlandford and Ostriker1ll978l ) on the accre- 
tion shock on the pole. 

Diffusive shock acceleration mechanism occurs 
on the collisionless-shock case, which will be real- 
ized when the kinetic energy of the plasma parti- 
cles exceeds the Coulomb potential of the plasma. 
The ratio between them are called as the binding 
parameter, (, and the collisionless shock occurs 
on the condition of £ < 1. The parameter can be 
described by 



C ~ 0.18 



10 15 cm- 



i/3 fkT plx 



V leV 



(7) 



where kT pTe is a temperature of the pre-shock mat- 
ters. Therefore, the accretion shock on the WD 
pole has a condition of a collisionless shock. In the 
po st-shock plasm a, according to equation (5.31) 
bv lSpitzer ( 1962 ) the collisional timescale between 



the electrons and ions. t; , is given as 



0.01 



kT sh \ 3/2 
17.6 keV ) 



(8) 

which is one order of magnitude smaller than the 
cooling time scale (typically ~ 0.3 sec) by equation 
(|4]). Therefore, the collisionless region exists at 
the top of the accretion column, just beneath the 
shock, although electrons themselves or protons 
themselves may be in the collisional equilibrium 



under high densities. This idea is supported by 
the fact that the shallower spin modulation was 
observed in the harder energy band (section 13. ip . 

We can assume that the magnetic field is par- 
allel to the shock normal. In this case, the accel- 
eration time-scale on the shock (t acc ) is derived 
as 

- 20 e E 
3 eBu s z 

where £, E, e, B are fluctuation of the magnetic 
field (> 1), the maximum energy of accelerated 
electrons, electric char ge, and the magnetic field 



stren g th, respectively ( Jokipiilll987[ iBamba et al 
120031 lYamazaki et al.l I2OO6I) Assuming £ 



1 

(Bohm limit), B = 20 MG, and u s at equation 
(H}, the t acc becomes 



= 2.3 x 10" 



—) 

1 TeVy 

-2 



B 



20 MG 



1.2 x 10° cm s" 



(10) 



which is very small compared with acceleration 
time scale in you ng supernova remnants (e.g., 



Bamba et al.l 120051 ). In the actual case of the ac- 



cretion flow onto the WD, the magnetic fluctua- 
tion may be £ 3> 1, but if it satisfies £ < 10 4 , the 
time scale of the acceleration t acc is much shorter 
than the time scale of the accretion i coo i ~ 0.3 s. 

The synchrotron loss time scale (i S ync) is, on 
the other hand, 



9.8 x 10 



-13 



E 



1 TeV 



B 



20 MG 



(ID 

from Rybicki & Lightman (1979). The maximum 



energy of accelerated electron is achieved when the 
acceleration time-scale becomes same to the syn- 
chrotron loss time-scale. The maximum energy of 
the accelerated electrons is, thus, 



E 



1 TeV 



2.1 x 10" 



-1/2 



B 



20 MG 



1.2 x 10 9 cm s" 



-1/2 



(12) 



Our result imply that the accretion shock of CVs 
can accelerate electrons up to ^GeV when £ = 
1, which is seen in the cases of young supernova 
remnants and pulsar wind nebulae. In cases of 
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£ 3> 1, which work out when the shock is close to 
the WD surface and the magnetic field dominates 
the dynamics of the flow, the maximum energy will 
be reduced as the dependency of £ -0 ' 5 ; i.e., E will 
be ~MeV at £ = 1000. Note that the maximum 
energy of protons can be ~ 1800 times larger than 
electrons, ~ GeV to TeV. The cut-off frequency of 
synchrotron X-rays (v cu t) is 

u cut = l.0xW-A_W-^J) (i\ 

cut V 1TeV / v 20MG / v 1 / 

= 27keV , 

which is higher than our observation band, under 
an assumption of the Bohm limit condition. 

5. Conclusion 

We have observed the polar, AM Her, in a very 
low state with Suzaku at the X-ray luminosity of 

L7 x 1029 (gr^) er s sec-1 in the °- 5 - 10 keV 

band. 

1. The object shows a flare phenomenon reach- 
ing the X-ray luminosity of 6.0 x 10 29 ( 91 d pc ) 

erg sec -1 in the 0.5 - 10 keV band, with 
a time scale of ~ 3700 sec (section 13.11) . 
The X-ray spectra can be reproduced by 
the thermal-plasma MEKAL model with the 
temperature of 8.67^^4 keV and the abun- 
dance of 0.76±g:|| solars (section 13.21) . In 
discussion on the time scale of the flare, we 
concluded that it should not be the cooling 
time scale of the plasma from constraints on 
the plasma size and density (section [4. II) . 

2. The X-ray light curves during the qui- 
escence show a clear spin modulation at 
0.1289273(2) days at BJD 2454771.581. 
These values are consistent with those in 
a high state (section l3.1|) . The spectra dur- 
ing the quiescence are well represented by 
MEKAL + PL model or a single CEMEKL 
model (section I3.3|) . The phase- resolved 
analyses also support these two models (sec- 
tion S3). 

3. From historical X-ray measurements of AM 
Her in various states, we found that the tem- 
perature is positively correlated with the vol- 
ume emission measure. The Suzaku results 



in a very low state also follow this relation. 
We made a rough estimation with a simple 
numerical calculation, and found that the 
fcT-VEM trend can be interpreted by a tra- 
jectory with a constant plasma shape in var- 
ious electron densities, as an one of possible 
solutions (section l4.2j) . 

4. The origin of the non-thermal emission was 
discussed, and a shock acceleration process 
at the shock front on the top of the accretion 
column was checked and proposed (section 
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